Electron Emission from Diamondoids: A Diffusion Quantum Monte Carlo Study 
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We present density-functional theory (DFT) and quantum Monte Carlo (QMC) calculations 
designed to resolve experimental and theoretical controversies over the optical properties of H- 
terminated C nanoparticles (diamondoids). The QMC results follow the trends of well-converged 
plane-wave DFT calculations for the size dependence of the optical gap, but they predict gaps that 
are 1-2 eV higher. They confirm that quantum confinement effects disappear in diamondoids larger 
than 1 nm, which have gaps below that of bulk diamond. Our QMC calculations predict a small 
exciton binding energy and a negative electron affinity (NEA) for diamondoids up to 1 nm, resulting 
from the delocalized nature of the lowest unoccupied molecular orbital. The NEA suggests a range 
of possible applications of diamondoids as low-voltage electron emitters. 



Designing the nanoscale building blocks for future nan- 
otechnological devices is one of the most rapidly evolving 
fields of materials science. In the area of nanoscale op- 
toelectronic materials, attention was originally focused 
on semiconductor nanoparticles constructed from Si [l| 
or Ge [3], as these nanoparticles can easily be inte- 
grated with existing Si fabrication techniques. Re- 
cently, however, it was discovered that large quanti- 
ties of high-quality, H-terminated C nanoparticles can 
be isolated from petroleum and separated into mono- 
disperse samples. j3| These diamondoids are expected to 
have several technologically useful optoelectronic prop- 
erties. In particular, the effect of quantum confine- 
ment is expected to push diamondoid optical gaps into 
the UV range, enabling a unique set of sensing applica- 
tions. Additionally, it has been demonstrated [J, [SJ that 
H-terminated diamond surfaces exhibit negative electron 
affinities (NEA's), suggesting that diamondoids will also 
have NEA's. This opens up the possibility of coating sur- 
faces with diamondoids to produce new electron-emission 
devices. 

The breakthrough in isolatiiig diamondoids has 
sparked interest in measuring 0, Q and calculating 
their structural and optoelectronic properties 0, M, [13 
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12l |. However, such studies have proved to be chal- 
lenging, and have produced several controversial results. 
In 1999, the results of X-ray absorption near-edge struc- 
ture (XANES) studies of diamond films produced by 
chemical vapor deposition (CVD) 6,] were used to infer 
the evolution of the nanoparticle gap with size. This 
showed a persistence of quantum confinement effects up 
to 27 nm, much longer than in Si or Ge nanoparticles, 
where quantum confinement effects disappear above 5- 
7nm. In contrast, recent near-edge absorption fine struc- 
ture (NEXAFS) studies of diamondoids prepared by hot- 
filament CVD[13] and high-explosive detonation wavesQ 
showed that quantum confinement effects disappear in 
particles larger than 4 nm. 

Theoretical models for the optoelectronic properties of 
diamondoids have also produced contradictory results. 
Raty et al. f?! studied the size dependence of the opti- 



cal gap using both time-independent and time-dependent 
density-functional theory (DFT) calculations, and con- 
cluded that quantum confinement effects disappear in 
nanoparticles larger than 1 nm. Raty et al. 0] also pre- 
dicted that the gaps of diamondoids with sizes between 
1 and 1.5 nm are below the gap of bulk diamond. This 
is strikingly different to the behavior of H-terminated 
Si and Ge nanoparticles, for which the gaps are always 
above the bulk band gap. 1J| In contrast, recent DFT cal- 
culations by Mcintosh et al. for the same particles 
predict optical gaps 2 eV above the gap of bulk diamond 
for particles ranging in size from 0.5 to 2 nm. 

In this Letter we resolve these controversies by per- 
forming both DFT and highly accurate quantum Monte 
Carlo (QMC) calculations to predict the size depen- 
dence of the optical gap and the electron affinity of di- 
amondoids. We have studied two classes of nanopar- 
ticles: (i) diamondoids constructed from adamantane 
cages, adamantane, CioHig, diamantane, C14H20, and 
pentamantane, C26H32; and (ii) H-terminated, spherical, 
diamond-structure nanoparticles, C29H36, C66II64, and 
C87H76 . The nanoparticle sizes and symmetries are given 
in Table HI Because diamondoids can be extracted in 
large quantities from petroleum and highly purified us- 
ing high-pressure liquid chromatography, we anticipate 
that actual experimental samples will consist largely of 
the high-symmetry structures studied in this work. This 
is not the case for Si and Ge nanoparticles, where limita- 
tions in current synthesis techniques prevent the routine 
production of high-symmetry nanoparticles. 

Our DFT calculations were performed using both 
plane-wave and Gaussian basis sets. The PBE [l6| func- 
tional was used in all the DFT calculations, as this has 
been shown to reproduce accurately the structural and 
electronic properties of diamond. [3] The pla ne- wave cal- 
culations were performed with the GP |l7l |. QBox[l8l|. 
and ABINIT [l9| codes. A plane-wave cutoff of 55 Ry 
and a cubic simulation cell of side-length 30 A were used 
for all calculations. Norm-conserving TrouUier-Martins 
(TM) pseudopotentials were used to represent the ionic 
cores. The all-electron Gaussian calculations were per- 
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Cluster 


Sym 


Size DFT gap 


DMC gap 


EA IP 


C10H16 


Ta 


0.50 


5.77 


7.61(2) 


-0.13(2) 10.15(3) 


C14H20 


Dm 


0.69 


5.41 


7.32(3) 




C26H32 


Ta 


0.74 


5.03 


7.04(6) 




C29H36 


Td 


0.76 


4.90 


6.67(6) 


-0.29(6) 7.63(5) 


C66H64 


Td 


1.00 


4.37 


5.09(18) 




C87H76 


Td 


1.14 


3.94 






Diamond 






4.23f7] 


5.6(2)fl5] 





Number of Carbon Atoms 
14 26 29 66 



TABLE I: Point-group symmetry, size (in nm), DFT and 
DMC optical gaps, DMC electron affinities (EA), and DMC 
ionization potentials (IP) (all in eV) of different diamondoids. 



formed with the G03 code.[20| 

The initial structures of CioHig, C14H20, and C26H32 
were taken from Ref. while those of C29H36, C66H64, 
and C87H76 were constructed from the bulk-diamond 
structure, and the C-H bond lengths were set to the ex- 
perimental value in CH4. All the structures were fully 
relaxed to their lowest-energy configurations within DFT 
before their electronic and optical properties were calcu- 
lated. For C10H16, C14H20, and C26H32, DFT-PBE cal- 
culations predict small relaxations in the bond lengths 
and angles, similar to Ref. [12] The C-C and C-H bond 
lengths are 1.55 and 1.10 A, respectively, essentially iden- 
tical to the bulk-diamond and methane bond lengths. 
The bond angles vary from 108.5° to 110.5°. For the 
larger, spherical particles (C29H36, C66H64, and CsrHre), 
DFT-PBE calculations predict a 3-4% increase in the C- 
C bond lengths, while the C-H bonds remain at 1.10 A. 
The bond angles vary from 102° to 113°. 

The results of our DFT calculations for the single- 
particle states of the diamondoids are shown in Fig. [T] and 
listed in Table [J Figure [T^ shows the size dependence of 
the DFT eigenvalue gap between the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO), calculated using the plane- 
wave basis set described above. All the HOMO-LUMO 
gaps are converged with respect to the size of the super- 
cell and the plane-wave basis. The HOMO-LUMO gap 
decreases as the size of the nanoparticles increases. For 
nanoparticles larger than 1 nm the DFT-PBE HOMO- 
LUMO gap is smaller than the DFT-PBE gap of bulk 
diamond. 

Figure [T]d shows the size dependence of the individ- 
ual DFT HOMO and LUMO eigenvalues. When the size 
of the nanoparticles is reduced, the HOMO eigenvalue 
decreases continuously, as predicted by standard mod- 
els of quantum confinement. 21| The lower symmetry of 
C14H20 does not alter the trend in the HOMO eigen- 
value. In contrast, the LUMO eigenvalue displays al- 
most no quantum confinement and is nearly independent 
of size. This behavior differs from that of H-terminated 
Si and Ge nanoparticles (l3| . in which both the HOMO 
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FIG. 1: Calculated size dependence of (a) DFT HOMO- 
LUMO gap and DMC optical gap, and (b) DFT HOMO and 
LUMO eigenvalues. The inset shows the convergence of the 
DFT HOMO and LUMO eigenvalues with basis set. 



and LUMO exhibit clear quantum confinement. 

The origin of the anomalous size dependence of the 
LUMO energy can be understood by comparing the 
HOMO and the LUMO. Figure [2] shows isosurfaces of the 
HOMO and LUMO of C29H36. The HOMO is locahzed 
on the C-C and C-H covalent bonds inside the nanopar- 
ticle, while the LUMO is a delocalized state with con- 
siderable charge outside of the H atoms terminating the 
surface. As the size of the nanoparticles increases, one 
expects the HOMO to evolve smoothly into the valence- 
band maximum of bulk diamond. In contrast, for larger 
diamondoids, the LUMO remains localized on the surface 
and is closer in nature to a defect or surface level within 
the gap of the bulk material; the LUMO does not evolve 
into the conduction-band minimum. This surface nature 
causes the optical gaps of the larger diamondoids to lie 
below the bulk gap. 

The surface nature of the LUMO is also responsible 
for some of the previous controversies in evaluating the 
energy of this state. The inset to Fig.[T)D shows an exam- 
ple of the convergence of the electronic eigenvalues in a 
DFT-PBE calculation of CioHig using a Gaussian basis. 
It shows the energy of the HOMO and LUMO calculated 
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FIG. 2: (Color) Isosurface plots of the square of the (a) 
HOMO and (b) LUMO of C29H36- The green isosurfaces in- 
clude 50% of the charge in each orbital. 

with the G03 code ^ using6-31G*, 6-311G**, aug- 
ccp-vtz,[22l and aug-ccp-vqz [2§| basis sets. The energy 
of the LUMO decreases by almost 2eV as more diffuse 
basis functions are added, while the HOMO energy stays 
approximately constant. If a localized basis set with in- 
sufhcient flexibility to describe the diffuse character is 
used, the LUMO will be artificially localized close to the 
nanoparticle, increasing its kinetic energy and pushing up 
its energy. However, when the large aug-ccp-vtz and 
aug-ccp-vqz [22,] basis sets are used, the Gaussian DFT- 
PBE HOMO-LUMO gaps agree with the plane-wave re- 
sults to within 0.01 eV. 

The QMC calculations were performed with the 
CASINO [i^l code using Slater-Jastrow trial wave 
functions of the form = exp[J], where 

and are Slater determinants of up- and down-spin 
orbitals taken from DFT calculations and exp[J] is a Jas- 
trow correlation factor, which includes electron-electron 
and electron-ion terms expanded in Chebyshev polyno- 
mials. The Jastrow factor was optimized using a stan- 
dard variance-minimization scheme, [i^l All the QMC op- 
tical ga ps were calculated using diffusion Monte Carlo 
(DMC).[2J| To converge the DMC total energies with re- 
spect to time step and population size, a time step of 
0.02 a. u. was used and the target population was at least 
640 configurations in each calculation. The ionic cores 
were represented by the same TM PBE pseudopotentials 
used in the DFT calculations. The effect on the QMC 
energies of using LDA or PBE functionals to generate the 
DFT orbitals was also tested. When the PBE fimctional 
and pseudopotentials were replaced with the LDA func- 
tional and pseudopotentials, the DFT and DMC gaps 
of C29H36 were reduced by 0.1 eV and 0.3(1) eV respec- 
tively. The DFT and DMC optical gaps are therefore 
relatively insensitive to the choice of functional. 

The optical gaps of the diamondoids were calculated 
as the difference in the DMC energy of the ground state 
and an excited state. The absorption of a photon cre- 
ates an excited singlet state. The description of such a 
state requires two pairs of Slater determinants. However, 



for computational simplicity, we represented the excited 
state by replacing the HOMO in the spin-down Slater de- 
terminant with the LUMO. Withiri_Ziegler's sum model 
for restricted Hartree-Fock states |25| , the error incurred 
by adopting this mixed-state approach is equal to half the 
singlet-triplet splitting. This error is typically 0. 1-0.2 eV 
in group- IV nanostructures of this size |2q | , which is small 
compared with the optical gaps. Previous QMC calcula- 
tions of the optical gaps of Si nanoparticles [IJ] using this 
approach were shown to be in excellent agreement with 
GM^-Bethe-Salpeter-equation (GW^-BSE) calculations of 
the true singlet excitation energies. 

To investigate the sensitivity of the QMC excita- 
tion energies to the choice of single-particle states, 
the excited-state calculations were repeated with three 
choices for the LUMO: (i) the unoccupied LUMO from 
a ground-state DFT calculation; (ii) the HOMO* occu- 
pied by the excited electron in a DFT calculation of an 
excited triplet; (iii) the HOMO* occupied by the excited 
electron in a DFT calculation of an excited mixed state. 
The DMC excited-state energies evaluated using (i) and 
(ii) were 0.3(1) eV lower than the energy evaluated us- 
ing (iii), and therefore we conclude that the excitation 
energy is relatively insensitive to the choice of DFT or- 
bital used to represent the excited electron in the Slater 
determinant. 

The DMC optical gaps of the diamondoids are given 
in Table HI and are plotted against nanoparticle size in 
Fig. [TJ DMC calculations do not suffer from the well- 
known DFT "band-gap problem," as they fully describe 
the interaction of the valence electrons with the elec- 
tron excited into the LUMO by the absorption of a 
photon, so that electron-hole correlation is accounted 
for. The DMC optical gaps are significantly larger 
than the DFT HOMO-LUMO gaps, as was found for Si 
nanoparticles. [3| For example, the DMC optical gap of 
C10H16 is 7.61(2) eV, while the DFT-PBE gap is 5.77 eV. 
Nevertheless, the DMC results confirm the qualitative 
trend for the size dependence of the gap predicted by 
DFT-PBE. The DMC calculations predict that quantum 
confinement will only be observed in diamondoids smaller 
than 1 nm. The rapid decay of quantum confinement in 
diamond nanoparticles is consistent with the small ex- 
citon Bohr radius in diamond (1.6 nm), compared to Si 
(4.9 nm). The optical gap of C66H64 is below the bulk- 
diamond gap, as Raty et al. Q predicted in their DFT 
study. These DMC optical-gap calculations support the 
interpretation of the X-ray absorption measurements in 
Ref. [rl, which found no quantum confinement effects in 
4nm diamondoids. 

To calculate the electron affinity {EA) and ioniza- 
tion potential {IP) of representative nanoparticles, DFT 
and QMC calculations were performed for the total en- 
ergy differences EA — En — E^+i and IP = Ejq-i — 
En, where E^ is the ground-state energy of a neutral 
molecule with N electrons. The trial wave functions for 
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the Em+1 and E^-i systems were constructed by adding 
and removing the orbitals corresponding to the LUMO 
and the HOMO from the down-spin Slater determinant. 

The DMC electron affinities and ionization energies of 
C10H16 and C29H36 are given in Tabic [D Both nanopar- 
ticles have a negative electron affinity (—0.13(2) and 
— 0.29(6) eV respectively). In contrast to the optical 
gaps, where the DMC values are significantly larger than 
the DFT values, the DMC electron affinities agree well 
with the DFT values; for example, the DFT electron 
affinity of CioHig is -0.128eV. The NEA suggests that 
coating surfaces with diamondoids could be a simple and 
economic method of producing electron emitters. 

Finally, we have studied the exciton binding energy, 
which is defined as the difference between the quasi- 
particle and optical gaps, for small diamondoids. For 
Si nanoparticles, QMC and GW-BSE calculations have 
found that the exciton binding is dramatically enhanced 
at the nanoscale due to the increased overlap of the elec- 
tron and hole. [13] Surprisingly, the exciton binding en- 
ergy of C29H36 is only 1.25(10) eV compared to 2.6eV 
for the equivalent Si29H3g nanoparticle. 27[ The origin of 
this reduced exciton binding can again be traced to the 
different nature of the LUMO in C and Si nanoparti- 
cles. The delocalized LUMO in C29H36 (see Fig.Eb) has 
a much smaller dipole matrix element with the HOMO 
than the equivalent orbital in Si29H3g, reducing the ex- 
citon binding. 

In conclusion, we have performed DFT and QMC cal- 
culations of the single-particle and optical gaps of H- 
terminated C nanoparticles. Our DFT and QMC cal- 
culations confirm the predictions of Raty et al. Q that 
quantum confinement effects disappear in nanoparticles 
larger than 1 nm. We have shown that the LUMO of H- 
terminated C nanoparticles is a delocalized surface state, 
in contrast to Si and Ge nanoparticles, where the LUMO 
is core-confined. The delocalized nature of the LUMO re- 
sults in an anomalously small exciton binding, a negative 
electron affinity, and optical gaps of larger nanoparticles 
that are below the bulk gap. 

We thank J. Grossman for performing the Gaussian 
DFT calculations, and J.-Y. Raty, T. van Buurcn, and 
J. Dahl for helpful conversations. N.D.D. acknowledges 
a fellowship in the LLNL Summer Institute on Computa- 
tional Chemistry and Materials Science. This work was 
performed under the auspices of the US Department of 



Energy by the University of California, Lawrence Liver- 
more National Laboratory under contract No. W-7405- 
Eng-48. N.D.D. and R.J.N, acknowledge financial sup- 
port from EPSRC, UK. Computational resources were 
provided by LLNL and the Cambridge-Cranficld HPCF. 



[2] 
[3] 

[5 

[6 
[7 

[s: 
[9: 

[10: 
[11 

[12 

[13 
[14 

[15 

[16 

[17 
[18 
[19 



[20 

[21 
[22' 
[23 

[24 
[23 

[26 

[27 



1242 



Electronic address: 'williamsonlO@llnl.gov' 
W. Wilson, P. Szajowski, and L. Brus, Science 262, 
(1993). 

G. Medeiros-Ribeiro et al, Science 279, 353 (1998). 

J. E. Dahl, S. G. Liu, and R. M. K. Carlson, Science 299, 

96 (2003). 

F. Himpsel et al, Phys. Rev. B 20, 624 (1979). 

M. J. Rutter and J. Robertson, Phys. Rev. B 57, 9241 
(1998). 

Y. K. Chang et al, Phys. Rev. Lett. 82, 5377 (1999). 
J.-Y. Raty et al, Phys. Rev. Lett. 90, 037401-1 (2003). 
S. Hosokawa et al., J. Elec. Spect. 137, 235 (2004). 
A. Barnard, S. Russo, and L Snook, Phys. Rev. B 68, 
073406 (2003). 

G. Lee et al, Phys. Rev. Lett. 91, 265701 (2003). 
N. Park et al, Phys. Rev. B 69, 195411 (2004). 

G. Mcintosh, M. Yoon, S. Berber, and D. Tomanek, 

Phys. Rev. B 70, 045401 (2004). 

Y. Tang et al, Chem. Phys. Lett. 372, 320 (2003). 

A. J. WiUiamson et al, Phys. Rev. Lett. 89, 196803-1 

(2002). 

M. D. Towler, R. Q. Hood, and R. J. Needs, Phys. Rev. 
B 62, 2330 (2000). 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 
F. Gygi, The GP Code, LLNL (2004). 
F. Gygi, The QBox Code, LLNL (2004). 
The ABINIT code is a common project of the Universite 
Catholique de Louvain, Corning Inc., and other contrib- 
utors, (http://www.abinit.org) (2004). 
M. Frisch, Gaussian 03, Revision A.l, Gaussian Inc., 
Pittsburgh, PA (2003). 

A. Yoffe, Advances in Physics 51, 799 (2002). 
T. J. Dunning, J. Chem. Phys. 90, 1007 (1989). 
R. J. Needs et al, CASINO version 1.6 User Manual, Uni- 
versity of Cambridge, Cambridge (2002). 
W. M. C. Foulkes et al. Rev. Mod. Phys. 73, 33 (2001). 
T. Ziegler, A. Rauk, and E. Baerends, J. Theoret. Chim. 
Acta 43, 261 (1977). 

D. Prendergast et al, J. Am. Chem. Soc. 126, 13827 
(2004). 

L. Benedict et al, Phys. Rev. B 68, 085310 (2003). 



